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Аннотация 

Введение. При бурении нефтяных или газовых скважин порода из забоя выносится на поверхность буровым рас- 
твором, который очищается от шлама на первой стадии виброситами. Вибросита оснащаются такой сеткой и 
создают такие траектории виброколебаний рамы, чтобы с одной стороны — обеспечивать очистку поступающего 
из устья скважины бурового раствора, а с другой — обеспечивать должную степень очистки. Для корректного 
подбора номенклатуры сеток нужно надежно определять пропускную способность вибросита в реальных усло- 
виях. Это позволит заранее заказать требуемое количество сеток с нужным размером ячеек. Ранее исследования 
проводились путем просеивания смеси калиброванных гранул на стандартном вибростенде или процеживанием 
минерального масла через неподвижную сетку. Однако это не в полной мере соответствует реальным условиям 
просеивания бурового раствора. Целью данной работы было моделирование пропускной способности бурового 
вибросита в реальных условиях путем расчета движения бурового раствора по вибрирующей сетке с использо- 
ванием её ранее полученной удельной пропускной способности. 

Материалы и методы. При создании математической модели полной пропускной способности бурового вибро- 
сита использовалась известная ранее математическая модель течения жидкости в открытом канале и конечно- 
мерная аппроксимация в виде последовательности участков из сосредоточенных емкостей, соединенных гидрав- 
лическими проводимостями (по аналогии с аппроксимацией электрических длинных линий и протяженных газо- 
проводов методом конечных элементов). Расход течения жидкости над сеткой определялся по закону Шези. При 
этом площадь сечения потока над сеткой для конкретного участка определялась как среднее арифметическое 
между начальным и конечным значением высоты бурового раствора на соответствующем участке. 

Результаты исследования. Составлена цепная расчетная схема просеивания бурового раствора. На основании 
материального баланса разработаны математические модели определения пропускной способности и высоты бу- 
рового раствора над сеткой: 

— для первого участка сетки и первой узловой точки сетки; 

— для последующих участков сетки и последующих узловых точек сетки; 

— для последнего участка сетки и крайней узловой точки сетки. 

Получена математическая модель движения бурового раствора по сетке вибросита, которая вследствие конечно- 
мерного представления этого потока по и участкам длины сетки представляет собой систему из и интегральных 
уравнений, описывающих сосредоточенные емкости, и п — | алгебраических уравнений, описывающие перетоки 
бурового раствора между емкостями. 

Обсуждение и заключение. Движение бурового раствора по длине сетки вибросита соответствует движению 
установившегося неравномерного потока в открытом русле. Однако в первом случае — объемный расход по 
длине переменный (убывающий до нуля), а во втором случае — постоянный по длине. Поэтому использование 
только одного дифференциального уравнения установившегося неравномерного потока в открытом русле в дан- 
ном случае недостаточно. 
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Полученная математическая модель движения бурового раствора по сетке вибросита является математическим 
обоснованием ДлЯ разработки моделирующей программы по определению пропускной способности вибросита 
для реальных условий его работы, то есть для конкретного бурового раствора, данной сетки или заданной си- 
стемы приводов вибросита. 
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Введение. Пропускная способность комплекса очистки бурового раствора от шлама циркуляционной си- 
стемы буровой установки лимитируется расходом циркулирующего бурового раствора, который должен быть 
очищен должным образом [1]. Под пропускной способностью понимается объем бурового раствора, просеянный 
сквозь сетку, а в качестве расхода — объем бурового раствора, протекший без просеивания. Вибросито является 
первым технологическим аппаратом в этом комплексе и осуществляет грубую очистку. Поэтому его пропускная 
способность в основном определяет пропускную способность всей системы очистки. 

Пропускная способность комплекса очистки циркуляционной системы (ЦС) стандартной буровой установки 
должна превышать производительность бурового насоса и составлять порядка 30-80 л/с. Недостаточная пропуск- 
ная способность ЦС ограничивает скорость проходки [2]. 

Циркулирующий буровой раствор поступает на вибросито напрямую из устья скважины по желобу. Здесь нет 
промежуточных емкостей, поэтому даже если его расход кратковременно превысит пропускную способность 
вибросита, то буровой раствор по сетке перельется в отвал, что приведет к его потере. Если помощник буриль- 
щика откроет шиберную заслонку, то неочищенный буровой раствор мимо сетки вибросита поступит в емкость. 
Это крайне нежелательный случай, так как шлам оседает в емкости, не позволяя перекачивать его насосом в 
пескоотделитель. Это резко снижает степень очистки бурового раствора, что приводит к увеличению затрат на 
бурение и даже к осложнениям в скважине [2]. 

Пропускная способность конкретного вибросита во многом определяется размером ячейки сетки, которой 
оно оснащено. С уменьшением размера ячейки существенно уменьшается пропускная способность вибросита, 
но повышается степень очистки бурового раствора. На практике для предотвращения пропуска раствора мимо 
вибросита на нем устанавливают сетку с более крупной ячейкой с большим, но неизвестным запасом по про- 
пускной способности [1]. 

Реальное измерение пропускной способности вибросита затруднено, т.к. после просеивания через сетку рас- 
твор сливается в емкость под виброситом самотеком, не проходя по трубопроводу. Поэтому определение про- 
пускной способности бурового вибросита проводилось косвенными методами. 

В работе [2] на основании производственного опыта указаны параметры, изменяющие пропускную способ- 
ность: увеличение длины сетки, снижение скорости потока, уменьшение угла наклона сетки, изменение направ- 
ления перемещения частиц, уменьшение амплитуды колебаний сетки, одновременное использование двух после- 
довательных или параллельных сеток. Однако авторы не приводят количественные показатели данных парамет- 
ров. В работе [3] приведены значения пропускной способности бурового вибросита при работе с различными 
буровыми растворами, но не приводится общая формула. В! указывается только формула для расчета изменения 
пропускной способности. Задача выбора сетки для конкретного случая решена в обзоре коллектива американских 
инженеров?. Она состоит в том, чтобы выбрать сетку с наименьшей ячейкой, при которой вибросито обладало 
бы пропускной способностью больше заданной. Для испытаний сеток разработано несколько методов, но все они 
далеки от условий работы вибросит на буровых установках. Так в стандарте АРТ Весопитепае4 ргасйсе 13С сетки 
испытываются или просеиванием смеси калиброванных гранул на стандартном вибростенде или процеживанием 
минерального масла через неподвижную сетку. В работе [4] сделана попытка на основе уравнения Навье-Стокса 
аналитически описать движение потока бурового раствора по сетке и процесс просеивания его через сетку. Ре- 
зультаты моделирования не учитывают динамические параметры грохота. Эти параметры учтены в работе [5], 
но в качестве исследуемого не для бурового раствора, а для угольного шлама. Количество всего удаляемого 
шлама от общей массы сравнительно невелико [6]. Такая концентрация шлама не оказывает заметного влияния 
на реологию бурового раствора [3]. Соответственно разделение процесса обезвоживания, приведенное в [7] не- 
целесообразно. Большая концентрация шлама наблюдается в буровом растворе, выходящего из пескоотделите- 
лей и илоотделителей, который очищается на втором вибросите?. 
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3 Там же. С. 3-335. 
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Из всего вышесказанного следует, что определить значение пропускной способности бурового вибросита с 
помощью теоретических исследований пока не получалось. Поэтому целью данной работы было моделирование 
полной пропускной способности бурового вибросита на основании проведенных ранее лабораторных экспери- 
ментов по определению удельной пропускной способности. 

Удельная пропускная способность. Коллективом с участием автора ранее была разработана методика [7], 
техническое и программное обеспечение для экспериментального определения математической модели удельной 
пропускной способности конкретной сетки первого вибросита по реальному буровому раствору для конкретного 
вибросита, проведены экспериментальные исследования [8]. Если у всех точек рамы сита, расположенных по 
линии сетки, траектории виброколебаний одинаковы, или, как говорят в вибротехнике «траектории однородны», 
то достаточно провести один эксперимент. Если траектории неоднородны, то достаточно провести три экспери- 
мента при траекториях на краях и в центре сетки. В любом случае параметры виброколебаний точек рамы у 
производственного вибросита должны быть известны и реализованы в экспериментах. 

Известны экспериментальные стенды других ученых, например, в работе [9] описан стенд, который создает 
полигармонические траектории колебаний. Это достигается соответствующей компоновкой системы привода с 
кратной частотной вращения. В ходе экспериментов получены даже фигуры Лиссажу [10]. Полученная наилуч- 
шая конструкция может работать менее эффективно за счет внешних воздействий. При этом целью автоматиче- 
ского управления является компенсация негативного влияния внешних воздействий. В работе [11] приведен спо- 
соб управления, в [12] — динамические параметры вибросита в ходе управления, а в [13] приведен пример по- 
шаговой оптимизации с помощью комплексов моделирования Ада и Мазгап. 

Проведенные эксперименты по просеиванию воднобентонитового раствора плотностью 1,19 г/см? через об- 
разец сетки с квадратными ячейками размером 0,4х0,4 мм позволили получить следующую математическую мо- 
дель удельной пропускной способности сетки [9] 


а(в, Ау) = 0,0088-й"8?. Ау", р. 


где 4(й, Ау) — удельная пропускная способность сетки, м?/м?с; й — высота слоя бутового раствора на сетке, м; 
Ау — амплитуда вертикального виброускорения, м/с?. 

Граничные условия процесса просеивания 

Наличие математической модели удельной пропускной способности недостаточно для определения пропуск- 
ной способности бурового вибросита. Необходимо также знание распределения высоты слоя бурового раствора 
по длине сетки й(х) при известном расходе раствора, поступающего на сетку сита 0(0). Принять в первом при- 
ближении линейное распределение высоты по длине сетки не представляется возможным, так как неизвестно 
граничное условие в конце процесса просеивания й(0) Граничное условие начала процесса просеивания, как было 
описано ранее, известно и имеет вид 1(0,8 . Г.) = 0, где Г — длина сетки, м. 

В классической гидравлике известна математическая модель распределения высоты слоя жидкости по длине 
потока при безнапорном установившемся неравномерном течении жидкости в открытом канале [14]. Однако рас- 
ход жидкости в этом случае постоянен по длине потока, а в нашем случае расход уменьшается вследствие про- 
сеивания части бурового раствора. Следовательно, установившееся движение потока бурового раствора по сетке 
вибросита представляет собой объект с распределенными по длине сетки переменными й(х) и О(х). 

Методы расчета распределенных объектов 

Поиск инженерных решений подобных задач показал, что в электротехнике электрические длинные линии 
заменяют цепочкой сосредоточенных сопротивлений, емкостей и индуктивностей [15]. Протяженные газопро- 
воды для моделирования их работы заменяют цепочкой перемежающихся емкостей и гидравлических сопротив- 
лений [16]. Особенно часто подобный подход используется в динамике машин, где по чертежам общих видов 
строятся цепные расчетные схемы, состоящие из перемежающихся сосредоточенных масс или моментов инерции 
и сосредоточенных податливостей для вибрационных машин [17] и систем управления [18]. Это по сути является 
конечно-элементной аппроксимацией распределенных объектов. 

Построение математической модели полной пропускной способности бурового вибросита 

Построим расчетную схему потока бурового раствора из сосредоточенных емкостей и сосредоточенных гид- 
равлических проводимостей, приведенную на рис. 1. Для этого разобъем длину сетки вибросита узловыми точ- 
ками на равные участки длиной Дх = [/т, где т — число участков. Узловые точки обязательно должны стоять в 
начале и в конце сетки. В узловые точки поместим сосредоточенные емкости. Начальной и конечной емкости 
соответствует участок потока длиной 0,5 . Ах, остальным емкостям — Дх. 
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О 1 0,1 О 1 


Рис. 1. Цепная расчетная схема гидравлической системы движения потока бурового раствора по сетке 


Уравнение материального баланса ДлЯ первой емкости на рис. 1 имеет вид: 


(2-0, ()-0. (1-4 Е -ь-ан (о), @) 


где О (В) — расход бурового раствора, поступающий на вибросито, м3/с; О5(1) — расход бурового раствора, 
вытекающий сетке из емкости 1-й узловой точки в емкость 2-й узловой точки, м3/с; О„1(В — расход бурового 
раствора, просеянный на первом участке, м3/с; й1(#) — высота слоя раствора в первой узловой точке, м; 
Ь — ширина сетки, м. 

Просеянный расход в первой узловой точке с учетом удельной пропускной способности равен: 


О (= ч(в (0.4). 5 (3) 


где Ау; — амплитуда виброускорения в первой узловой точке. 
Подставив выражение (3) в (2) и преобразовав уравнение, получим выражение для й1(1: 


2 
(д [оо - ода. @ 
Ах.Б 
Уравнение материального баланса для емкостей всех промежуточных узловых точек на рис. 1 имеет вид: 

(0. (1-0 (1-0 (1))-# = к-ь-4 (0, (5) 
где О{Е) — расход бурового раствора, втекающий в емкость {-той узловой точки, м3/с; О;--1(8) — расход бурового 
раствора, вытекающий по сетке из емкости {-й узловой точки в емкость /+1-й узловой точки, м3/с; 
О„Ка) = а(#ко, Ау) : Ах : Б — расход бурового раствора, просеянный на {-ом участке, м/с; й(#) — высота слоя рас- 


твора в 1-ой узловой точке, м. 
Преобразовав уравнение (5), получим выражение для 11) всех промежуточных узловых точек: 


1; (1) = = (0:(1)-О-(1))-а (№: (1), Аз )|-4 (6) 


В последней узловой точке исходя из физических соображений примем высоту слоя раствора нулевой (С), 
чтобы не подпирать поток в случае течения его в отвал. Это по сути — правое граничное условие для высоты 
слоя раствора. Расход бурового раствора, поступающий на вибросито, в данной математической модели принят 
постоянным и равным заданной величине 01(1) = О! зал. Это косвенное задание левого граничного условия для 
высоты слоя раствора. 

Все промежуточные расходы ОК® на рис. 1 определяются перепадом высот слоя раствора й: 1(1) — #1) и сосредото- 
ченным гидравлическим сопротивлением участка сетки. Буровой раствор является неньютоновской жидкостью, но 
при турбулентном течении глинистых растворов можно исходить из обычных формул гидравлики [19]. 

Для определения расхода Ой не будем определять величину гидравлического сопротивления участка сетки 
длиной Дх, а воспользуемся известным дифференциальным уравнением установившегося неравномерного по- 
тока постоянного расхода [20], но для поднятой сетки: 


аи (х) _ (+. 0? у. =. (1) 
ах Е?.С?.В в ? ] 


где 1. — геометрический уклон русла потока, на рисунке 1 й.> 01, >0; Е — площадь поперечного сечения по- 


тока, м?; С — коэффициент Шези; К — гидравлический радиус сечения потока, м; а — коэффициент Кориолиса, 
а = 1,10 - 1,15 [21]; ® — ускорение свободного падения. 


Кичкарь И.Ю. Математическая модель движения бурового раствора по сетке вибросита 


Выразим в дифференциальном уравнении (7) производную конечными разностями и подставим в него пере- 
менные из расчетной схемы на рис. 1: 


00 
Из (1) 1 (Е) “ Е2.С2.В, 
== 2 › (8) 
Дх о О; (1) 
о 
& В 


где Ё; — площадь поперечного сечения потока на участке от 1 - 1-ой узловой точки до 1-ой узловой точки. 

Здесь следует уточнить, что О0{) — это расход на участке от /— 1-ой узловой точки до 1-ой узловой точки. 
Площадь поперечного сечения потока на этом участке изменяется от величины й: 1(1) 6 до величины й® Ь. По- 
этому величину РА; будем определять как среднее значение этих площадей: 


(И (и). 


Е =5. 9 
2 (9) 
Величина гидравлического радиуса потока равна отношению [15]: 
т ЙЕ); (1 
т со 
7 Ь +2. 1-1 Й 


2 
где А; — длина смоченного периметра потока на участке от 1 — 1-ой узловой точки до 1-ой узловой точки, м. 
Ширина сетки промышленных вибросит — от 1 метра, толщина слоя раствора на сетке — несколько санти- 
метров [4]. Поэтому из выражения (10) следует, что А; = (й: (8) + #К0)/2. Подставив это значение и значение А’; из 
выражения (9), получим: 


+ ог (2) 
в (6)-в (6) (0,5-(в- (+ (2) -5?-с? в 
- р ог (1) | 


8 (0,5- (1-1 (19+ (1) -5° 
Разрешим уравнение (11) относительно ОКВ для {> 1: 
1 (1): (Е) ей 


2 Е г +в 5. () 


Ах ‘= с” 


Обсуждение и заключение. Математическая модель течения бурового раствора по сетке вибросита представ- 
ляет собой систему из и уравнений (6) и п- 1уравнений (12). 

Начальные условия для переменных й0) принять нулевыми нельзя, так как это приведет к делению на ноль 
в уравнениях (12). Поэтому рекомендуется принять эти условия ненулевыми, но достаточно малыми, например, 
#0) = 0,0001. 

Моделирование системы уравнений (6) и (12) можно производить в моделирующих средах, например, 
МАТГАВ с приложением ЗипайпК или ЗиишТесв. В любом случае надо определиться с количеством узловых 
точек цепной расчетной схемы на рис. 1. Пропускная способность вибросита, оснащенного конкретной сеткой, 
равна расходу О1 зал, при котором последние 20-25 % длины сетки были бы свободны от бурового раствора и 
использовались для обезвоживания шлама“. Поэтому предпоследняя узловая точка при минимальном количестве 
этих точек должна находится на таком же расстоянии от края сетки. Отсюда следует Ах = 0,15 : Г, чему соответ- 
ствует семь узловых точек. 

В этой работе стояла задача определения установившегося движения потока бурового раствора, а точнее 
его профиля. Переходный процесс этого движения автор вынужден был использовать, так как невозможно 
было установить левое граничное условие й1(с0), вместо которого использовано косвенное задание этого гра- 
ничного условия (01(0) = О! зад. Окончательное количество узловых точек определится только в процессе моде- 
лирования по показателю точности профиля бурового раствора. Для этого производят моделирование матема- 
тической модели (6) и (12) спил + 1 узловыми точками. Если наибольшая разность профилей потока превос- 
ходит, например, 5 % от значения высоты слоя раствора в этой точке, то увеличивают число узловых точек на 
единицу дои-2 [22]. Затем производят моделирование с и + 2-ми узловыми точками и сравнивают профили 


* Атенсап Аззостаноп о ОгИйтз Епятеегу ре ракету ап ОгИйиз Е 5уяет. Нопзюп: Ош Рите Сотрапу; 1999. 335 р. 
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потока этого моделирования с результатом предыдущего моделирования с п + 1 точками. По достижению за- 
данной погрешности в дальнейшем используют математическую модель с последним значением количества 
узловых точек. 

Полученная математическая модель движения потока бурового раствора по сетке вибросита совместно с ра- 
нее полученной математической моделью удельной пропускной способности сетки по конкретному буровому 
раствору [8] позволяет обоснованно рекомендовать установку сетки с определенными размерами ячеек для раз- 
личных расходов бурового раствора, истекающего из устья скважины. 
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